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ABSTRACT 

We present evidence for approximately 400-d variations in the radial velocity of HD 154791 (V934 
Her), the suggested optical counterpart of 4U 1700+24. The variations are correlated with the previously 
reported « 400 d variations in the X-ray flux of 4U 1700+24, which supports the association of these two 
objects, as well as the identification of this system as the second known X-ray binary in which a neutron 
star accretes from the wind of a red giant. The HD 154791 radial velocity variations can be fit with an 
eccentric orbit with period 404 + 3 d, amplitude K — 0.75 + 0.12 kms~^ and eccentricity e = 0.26 + 0.15. 
There are also indications of variations on longer time scales > 2000 d. We have re-examined all 
available ASM data following an unusually large X-ray outburst in 1997-98, and confirm that the 1-d 
averaged 2-10 keV X-ray fiux from 4U 1700+24 is modulated with a period of 400 + 20 d. The mean 
profile of the persistent X-ray variations was approximately sinusoidal, with an amplitude of 0.108 + 
0.012 ASM count s^^ (corresponding to 31% rms). The epoch of X-ray maximum was approximately 
40 d after the time of periastron according to the eccentric orbital fit. If the 400-d oscillations from 
HD 154791/4U 1700+24 are due to orbital motion, then the system parameters are probably close to 
those of the only other neutron-star symbiotic-like binary, GX 1+4. We discuss the similarities and 
differences between these two systems. 

Subject headings: binaries: symbiotic — X-rays: binaries — stars: neutron 



1. INTRODUCTION 

The X-ray source 4U 1700+24 (also known as 2A 
1704+241; I = 45?15, h = 32?99) was first discovered 
in Ariel V scans for high-latitude X-ray sources (Cooke 
et al. 1978), and by the Uhuru (SAS A) X-ray observa- 
tory (Forman et al. 1978). It has a typical flux of (1- 
10) X 10"" ergcm-^s"^ (2-10 keV, equivalent to 0.5- 
5 mCrab as measured by a range of X-ray missions; 
Masetti et al. 2002), and experiences occasional outbursts. 
During a 100-d X-ray high state in 1997-98, 4U 1700+24 
was observed by the Rossi X-ray Timing Explorer [RXTE: 
Jahoda et al. 1996). This high state was the larger of two 
known outbursts of this system, and the X-ray flux reached 
a peak of 40 mCrab. The X-ray spectrum of 4U 1700+24 
is similar to that of other accreting neutron stars, and 
generally requires a blackbody with a temperature of 0.9- 
1.3 keV plus a hard component at higher energies for an 
acceptable flt. The X-ray spectral hardness and the lack 
of UV continuum imply that the compact object is almost 
certainly not a white dwarf (Garcia et al. 1983). A black 
hole companion is also unlikely since the thermal X-ray 
component in black hole candidates is typically an order 
of magnitude lower in temperature than that measured in 
4U 1700+24. However, no rapid quasi-periodic or coher- 
ent X-ray variations have been confirmed (Masetti et al. 
2002), which is rather unusual for an accreting neutron 
star. 

Based on Einstein position measurements, Garcia et al. 
(1983) proposed the association of 4U 1700+24 with the 



8th magnitude M2 III star HD 154791. In a recent re- 
analysis of ROSAT HRI data, however, Morgan & Gar- 
cia (2001) found only a 10% chance that the two sources 
are associated, but did not propose an alternative candi- 
date. Previous optical spectroscopy ruled out radial ve- 
locity variations with amplitudes greater than 5 kms~^ 
on time scales of 40 min to ~ 1 yr (Garcia et al. 1983), 
strongly suggesting that either the orbital period is signif- 
icantly longer than one year, or that the system is viewed 
very close to face-on (inclination angle i w 0). 

HD 154791 may be only the second symbiotic-like binary 
known to contain a neutron star, after the well-studied 
M-giant/X-ray pulsar system V2116 Oph/GX 1+4. How- 
ever, these systems are quite different in X-rays (Table 1). 
Furthermore, the optical spectrum of V2116 Oph exhibits 
extremely strong, variable emission lines, probably pow- 
ered by UV photons originating from an accretion disk 
(Chakrabarty & Roche 1997). In contrast, HD 154791 has 
much weaker line emission. In fact, optically the spectrum 
is very close to that of an isolated M giant (e.g. Garcia 
et al. 1983). 

Here we describe optical and X-ray observations 
which go some way towards resolving the uncer- 
tainties surrounding this system. We compare the 
properties of HD 154791/4U 1700+24 with those of 
V2116 Oph/GX 1+4, and discuss what may be inferred 
about the former in light of our results. 
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Orbital phase 

Fig. 1. — Radial velocity measurements of HD 154791 and the 
corresponding 404 ± 3 d orbital solution. The top panel shows the 
radial velocity measurements versus time and the estimated la un- 
certainties. The bottom panel shows the measurements folded on 
the orbital period. In both panels the eccentric orbital model is 
overplotted as a solid line (see Table 2). 

2. OBSERVATIONS 

2.1. Radial velocity variations 

Echelle spectra of HD 154791 were obtained using inten- 
sified Reticon detectors on the 1.5-m Tillinghast Reflector 
at the Whipple Observatory on Mt. Hopkins, Arizona, 
and also the 1.5-m Wyeth reflector at Oak Ridge Observa- 
tory, Massachusetts. The observations were made between 
1982 April 12 and 1996 October 26, with the 18 observa- 
tions before 1992 October 7 taken at Mt. Hopkins, and 
those since taken at Oak Ridge. The instruments are de- 
scribed in detail by Latham (1992). The spectra cover a 
44 A bandpass centered around 5200 A, with a resolution 
of approximately 8.5 kms~^. Each observation consisted 
of a 10 20 minute integration of HD 154791, with 90 s 
Th-Ar comparison scans before and after. The positions 
of around 30 strong spectral lines in each calibration scan 
were fit with a fifth-order polynomial, and the pre- and 
post-observation scans were then combined to calibrate 
the wavelength scale for each spectrum. Typical residu- 
als to the polynomial fit were 0.5 pixel, corresponding to 
0.01 A, or 0.5 kms~^. Approximately 25 lines were com- 
mon to the entire set of wavelength solutions. The radial 
velocity was measured by cross-correlating each observed 
spectrum with an M-star template spectrum (e.g. Tonry 
& Davis 1979) and measuring the mean offset. Quadratic 
least-squares fits to the peak in the cross-correlation func- 
tion give a radial velocity accurate to better than 1 kms"-'^. 

The radial velocity measurements show significant vari- 
ations of approximately 1 kms~^ about the mean sys- 
temic velocity of —48.7 ± 0.1 kms~^ (Figure 1). We 
searched the data for periodic signals in the range 50- 
1000 d by calculating the Lomb-Scargle periodogram (e.g. 
Press et al. 1996), and found a marginally significant peak 
near a period of 410 d (Lomb-normalized power 10.5; esti- 
mated significance from bootstrap simulations, 3.3cr). Fit- 
ting an elliptical orbit to the data gave a solution with 
orbital period 404 ± 3 d and radial velocity amplitude 
K = 0.75 ± 0.12 kms-i (Table 2; solid lines in Figure 
1). The reduced-x^ (xt) ^ov the fit was 1.26, indicating a 



marginally acceptable fit. While the best-fit eccentricity 
was 0.26, the deviations from circularity were significant 
at less than the 2a level. 

There was also evidence for variations in the radial ve- 
locity on longer time scales. We calculated xt for an ec- 
centric orbit as a function of trial period for periods up to 
4000 d. While the 404 d signal gave rise to a xi value close 
to 1, lower values were achieved at substantially longer 
periods around 1930 d (xl = 1-04), 2370 d (1.02) and 
3700 d (0.91). The irregular spacing of the measurements 
did not allow us to distinguish clearly between these pe- 
riodicities; eccentric orbital solutions with any one of the 
three periods can account for the long-term radial veloc- 
ity variations. The observations do not cover more than 
one full cycle for these longer-period signals, so more ob- 
servations are needed to determine whether they are true 
periodic variations or aperiodic trends related to irregular 
red-giant variability. 

2.2. RXTE/ASM observations 

Masetti et al. (2002) noted tentative evidence for an 

approximately 400 d periodicity in RXTE All-Sky Mon- 
itor (ASM) data. To investigate the 400-d periodic- 
ity, we obtained 1-d averaged ASM measurements of 
4U 1700-h24 between 1996 January 5 and 2002 April 
18 (MJD 50087-52382) from the ASM WWW site at 
http://xte.mit.edu. The mean quiescent ASM count 
rate was 0.2 count s^^. During the flare observed between 
1997 September-December (MJD 50700-50800), the 1-d 
averaged rate reached 2.5 counts"^ (Figure 2, top panel). 

To avoid a biased period measurement due to the pres- 
ence of the 1997 98 outburst, we initially calculated a 
Lomb-normalized periodogram over the interval from the 
end of the outburst (around 1998 March 28) to 2002 April 
18. The result was a peak with Lomb-normalized power 
14.05 (equivalent to 3.9cr), at 383 ± 30 d. We estimated 
the period uncertainty by folding the ASM data on the 
detected period, and calculating 1000 sets of statistically 
equivalent data sets with the same time sampling, aver- 
age proflle and noise properties per bin as measured in 
the original data set. The variance in the detection pe- 
riods of the simulated data sets is then a measure of the 
uncertainty in the actual detection. 

To obtain a more precise estimate of the quiescent X- 
ray periodicity, we excluded those measurements with er- 
rors of > 0.5 count s""'^. Around 13% of the 1-d averaged 
ASM measurements had errors above this thrcsliokl, pos- 
sibly due to source variability on time scales less than a 
day, but more likely arising from low-significance detec- 
tions due to solar contamination or other sub-optimal ob- 
serving conditions. These low-significance measurements 
appeared to introduce a bias to the period measurements, 
and also tended to contribute disproportionately to the 
estimated uncertainty in the period. A Lomb-normalized 
periodogram calculated on this subset of data following 
the outburst resulted in a peak power of 24.04 (estimated 
significance 4.2 x 10~^, equivalent to 5.5(t) at 392 ± 14 d. 
To confirm the above results, we used the phase-dispersion 
minimization technique of Stellingwerf (1978) and folded 
the ASM light curve on a grid of periods between 200 to 
1000 d. We found that the X-ray fiux was most likely mod- 
ulated with a period of 404 ± 13 d. Combining the results 
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from the latter two period-search techniques, we determine 
the period of the X-ray oscillation to be 400 ± 20 d with 
the epoch of maximum at MJD 51140 ± 10 d. The mean 
profile of the post-outburst ASM measurements folded on 
the 400 d period is also shown in Figure 2 (inset, bottom 
panel). The profile was approximately sinusoidal with an 
amplitude of 0.108 ± 0.012 ASM count s-\ or 31% rms. 
With the periastron time from the orbital fit used as a ref- 
erence phase for the fold, we find that on average the X-ray 
maximum occurs 0.1 in phase 40 d) after periastron. 

For comparison, we also calculated an oversampled 
Lomb-normalized periodogram of the entire ASM data set 
(including the 1997-98 outburst), which revealed a sig- 
nificant periodicity at 383 d (identical to that measured 
following the outburst). The 1997-98 outburst was un- 
usual; no similar events were observed by the ASM in 
« 6.5 yr of monitoring. Thus, it may have occurred due 
to a mechanism distinct to some extent from the periodic 
behavior of the quiescent emission. Furthermore, there 
is evidence that the behaviour of the persistent emission 
prior to the outburst was different than afterwards. Al- 
though only a little more than one 400-d interval was cov- 
ered before the outburst, the persistent flux was weakly 
anticorrelated compared to the mean profile measured fol- 
lowing the outburst. Applying the same error threshold 
((7 < 0.5 counts"^) to the full data set resulted in an 
increase in both the peak power and period, to 51.1 (es- 



timated significance 1.2 x 10 



> 8cr) at 385 d. The 



shorter measured periods using the full ASM data set may 
be attributed to aliasing resulting from a delay between the 
nearest phase of X-ray maximum for the persistent emis- 
sion and the peak of the 1997-98 outburst. If 4U 1700-^24 
is similar to wind- accreting neutron stars in high-mass 
binaries, we may expect large outbursts to correspond 
only approximately with periastron phase (e.g. Stella et al. 
1986). In fact, the peak of the 1997 outburst falls around 
MJD 50780, approximately 40 d after the nearest epoch 
of maximum flux according to the 400-d periodicity de- 
tected in the post-1997 outburst ASM measurements, and 
approximately 70 days after periastron passage according 
to the orbital solution discussed in §2.1 (Figure 2). 

3. DISCUSSION 

We have found evidence for significant variation in the 
radial velocity measurements from HD 154791, which we 
fit with an eccentric orbit with period 404 ± 3 d. We 
have also found evidence for a significant variation in the 
X-ray flux from 4U 1700-1-24. Our best estimate of the 
X-ray period is 400 ± 20 d. Hence, the periods measured 
from the radial velocity and ASM X-ray flux variations 
are consistent. It is compelling to suggest that these two 
periodicities arise from a common source. Naturally, it 
is also possible that the agreement between the two peri- 
odicities is merely coincidence, and that HD 154791 and 
4U 1700+24 arc imrelated. Estimating the probability of 
such a coincidence depends upon the origin of the period- 
icities in each star, as well as the intrinsic distributions of 
these phenomena. In practice, it is not possible to make 
such an estimate with much confidence. However, it is 
clear that the likelihood of an unrelated red giant and X- 
ray source being at nearly the same measured position at 
high galactic latitude, and also both having a periodicity 
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Fig. 2.— RXTE/ ASM measurements of 4U 1700-1-24. The top 
panel shows the 1 d averaged ASM measurements between 1996 
January 5 and 2002 April 18 (MJD 50087-52382). The horizontal 
lines ( ) bracket the estimated times of phase maximum for the 

best period estimate (400 d; see §2.2). The open triangles indicate 
the predicted times of periastron passage for the eccentric orbital 
fit (^orb = 404 d; Table 2). The bottom panel shows the peri- 
odogram calculated from a subset of the ASM data after the end 
of the outburst (1998 March 28), where the measurement error was 
< 0.5 count s~^. The dashed vertical line shows the candidate 400 d 
period. The inset shows the folded ASM profile using the time of 
periastron from the orbital fit as the reference phase. 

that would be unusual for either type of source separately, 
is extremely small. The detection of a 400-d modulation 
from both stars therefore supports the original suggestion 
that HD 154791 is the optical counterpart of 4U 1700-1-24. 

3.1. Cause of the Correlation 

There arc several plausible mechanisms for a correspon- 
dence between the radial velocity and the persistent X-ray 
flux. We consider two possibilities for the origin of the 404- 
d period in the HD 154791 radial velocity curve: orbital 
motion, and pulsational variation of the M-giant. Precise 
optical photometric monitoring could distinguish between 
these two possibilities, since the maximum brightness oc- 
curring at any phase besides that of maximum recessional 
velocity would favour orbital motion. 

If the radial velocity variations are due to orbital motion, 
the 400-d X-ray flux modulation may arise from changes 
in the accretion rate onto the neutron star as it follows an 
elliptical orbit around the companion and moves through 
different nebular densities. Such modulation is common in 
high-mass wind- accreting X-ray binaries (e.g. Stella ct al. 
1986). Moreover, the times of periastron passage in our 
best-fit orbital model are very well correlated with the 
times of increased X-ray flux. For a 1.4 M© neutron star 
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and red giant mass < 2.5 Mq, the 404 d period and low 
amplitude of the radial velocity variations imply that the 
inclination is i < 2°. Although the a priori probability of 
observing such a system is < 7x lO""*, such a low incli- 
nation could explain the lack of X-ray pulsations or QPOs 
from the source (see §3.2). If the 400-d variation is in fact 
due to orbital motion, the longer-term variations in the 
radial velocities may be related to irregular variations in 
the red giant (Percy et al. 2001). 

Alternatively, red-giant pulsations with P = 400 d could 
modulate the stellar wind, and hence the accretion rate 
onto the neutron star. The IRAS fluxes of HD 154791 
(/i2 = 2.77, /25 = 0.736, feo < 0.4, /loo < 1.0, where 
/i2 is the flux density at 12 microns in Jy) indicate 
that is it not a Mira variable (Kcnyon et al. 1988). How- 
ever, non-Mira M-giants also pulsate. Although typical 
fundamental pulsation periods lie in the range 20-200 d, 
variations with a period up to an order of magnitude larger 
are also seen (Percy et al. 2001). If the 400-d variations 
are due to pulsation, one of the slower trends observed in 
the radial velocity variations could be related to orbital 
motion. In this case, the possible range of orbital incli- 
nations becomes much more likely. An orbital period of 
2370 d gives an inclination of 4-6° while for a period of 
3700 d the likely range is 5-7°. Also, the low accretion rate 
onto the neutron star would be more understandable if the 
orbital period is longer than 400 d. However, a longer or- 
bital period alone cannot account for the discrepancy in 
luminosity between 4U 1700-^24 and GX l-t-4. 

Finally, the shape of the 400-d radial velocity variation 
argues against a red-giant pulsation. Pulsation light curves 
are typically non-sinusoidal (in the sense that more time is 
spent with brightness below the median value than above 
it), and similar behavior would be expected for the radial 
velocity variations. 

3.2. Comparison with GX 1+4 

If the 400 d periodicity in the HD 154791 radial veloc- 
ities is due to orbital motion, then the orbital periods of 

HD 154791/4U 1700+24 and GX H-4 may be quite similar 
(see Table 1). This similarity makes the contrast between 
their other observed properties puzzling. The X-ray lu- 
minosity of GX 1-1-4 appears to be two to three orders of 
magnitude greater than that of 4U 1700-^24. In GX 1-1-4, 
there is clear evidence for an accretion disk (Chakrabarty 
& Roche 1997), whereas in 4U 1700-1-24 there is no mea- 
surable UV continuum (Garcia et al. 1983). GX 1+4 has a 
rich optical emission-line spectrum, and HD 154791 shows 
very little in the way of emission lines. Finally, GX 1+4 
is a 2-minutc X-ray pulsar, whereas no coherent or quasi- 
period oscillations have been convincingly detected from 
4U 1700+24. 

One explanation for the differing X-ray luminosities is a 
larger mass-loss rate from the late- type giant in GX 1+4. 
Chakrabarty & Roche (1997) find that the M6 giant in 
GX 1+4 is probably near the tip of the first-ascent red- 
giant branch. HD 154791 has a spectral type M2, and so 
would be expected to be losing mass at a lower rate. The 
relative lack of optical/UV emission lines suggests that 
either the ionized nebula in HD 154791/4U 1700+24 is 
substantially smaller and/or less dense than in GX 1+4, 
or there are not enough UV photons to power nebular line 



emission in 4U 1700+24/HD 154791. Either of these con- 
ditions can plausibly arise if the mass loss from the com- 
panion is much smaller. Moreover, if the binary separation 
is smaller in GX 1+4 than in 4U 1700+24/HD 154791 a 
larger fraction of the red giant wind could be accreted. 
The higher X-ray luminosity in GX 1+4 could also illumi- 
nate the red giant and increase the mass loss from the red 
giant further, in a sense producing a feedback effect. The 
observed variability of the optical spectrum of V2116 Oph 
is evidence for the significant role played by X-ray heating. 
On the other hand, in 4U 1700+24/HD 154791 the X-ray 
luminosity is perhaps low enough that very little illumi- 
nation of the red giant or enhanced mass loss is expected. 
Even during the 1997 X-ray outburst, Tomasella et al. 
(1997) note that the optical spectrum was unchanged. 

Finally, the lack of pulsations or QPOs in 4U 1700+24 
remains puzzling. Garcia et al. (1983) and Morgan & 
Garcia (2001) proposed a 900-s QPO, but this QPO was 
not confirmed in high-quality RXTE data (Masetti ct al. 
2002). High-precision optical photometry has also revealed 
a lack of pulsations (down to a limit of 0.5% fractional 
variation), or rapid variations of any kind (Sokoloski et al. 
2001). This behavior is again in contrast to GX 1+4, 
where Jablonski et al. (1997) find an oscillation in the op- 
tical emission with an amplitude of from 0.4 to 4.5%, plus 
stochastic variations. Given that the optical light from 
GX 1+4 has a significant contribution from an accretion 
disk, the lack of rapid optical variations from HD 154791 
provides evidence that any contribution from a disk in 
4U 1700+24 is negligible in the optical regime as well as 
the UV regime. The lack of X-ray pulsations is expected if 
the 400-d radial velocity variations arc due to orbital mo- 
tion. As discussed in §3.1, the low amplitude of the radial 
velocity variations requires i ~ 0, in which case any mag- 
netic hot spots on the neutron star will be continuously in 
view (assuming the NS spin axis is aligned with angular 
momentum axis of the binary). For higher and a priori 
more probable values of i, the non-detection of pulsations 
may be explained if the accretion flow is not appreciably 
disrupted by the magnetic field of the neutron star above 
the surface (e.g. Ghosh & Lamb 1979a, b), and spreads the 
accreting material more or less evenly over the neutron 
star surface. Given the extremely low luminosity (and 
hence accretion rate) of this source, standard calculations 
of ram and S-field pressure balance would require that 
the surface magnetic field strength of the neutron star be 
< 10^ G, far below the canonical value for low-mass X-ray 
binaries or even old neutron stars in recycled ms pulsars. 
A more likely scenario consistent with larger i is that the 
neutron star spin is very slow. 

In conclusion, all scenarios for 4U 1700+24/HD 154791 
have some difficulties. We have presented some options, 
but more observational constraints are needed to select be- 
ween them. Additional observations to determine whether 
the 400-d radial velocity variations are due to orbital mo- 
tion or red-giant pulsation would also enable one to better 
interpret the differences between 4U 1700+24/HD 154791 
and GX 1+4. 
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Table 1 

Comparison of the properties of the two known symbiotic neutron-star binaries. 



Parameter 


4U 1700+24 


GX 1+4 


References 


Distance (kpc) 


0.42 ±0.2 


3 6/12 15 


[1,2] 


-^spin (s) 


? 


« 100-137.7 


[3,4] 


Porh (d) 


404 ± 3? 


304? 


[5,6,7] 


Lx (2-10 keV, lO^^ ergs'^) 


(0.2-10) X 10-3 


~ 10/100'' 


[1,2] 




(0.1-5) X 10-3 


5/50*^ 




Companion spectra 


M2 III 


M3-6 III 


[1,2] 


Lopt/Lx 


200 


0.25 


[1] 



"^Depending upon the evolutionary status of the mass donor, first ascent red 
giant/beginning inital ascent of the asymptotic giant branch. 

''Depending upon the distance 

'^Assuming a neutron star with ii* = 10 km and M = 1.4M0. 



References. — 1. Masetti et al. (2002); 2. Chakrabarty & Roche (1997); 3. 
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Table 2 

Orbital parameters for HD 154791 

Radial velocity amplitude, K (kms-^) 0.75 ± 0.12 

Projected scmimajor axis, ax sin i (10^ km) 4.2 ± 0.7 

Orbital period, Porb (d) 404 ± 3 

Epoch of periastron, Tq (MJD) 49090 ± 80 

Eccentricity, e 0.26 ± 0.15 

Longitude of periastron, co (degrees) 260 ± 40 

Mass function, /o (10-^ Mq) 1.8 ± 0.9 

Model fit 95.96 (76 dof) 



